FINITE MASS SELF-SIMILAR BLOWING-UP SOLUTIONS OF A 
CHEMOTAXIS SYSTEM WITH NON-LINEAR DIFFUSION 



AD MEN BLANCHET 1 AND PHILIPPE LAURENgOT 2 



Abstract. For a specific choice of the diffusion, the parabolic-elliptic Patlak- 
Keller-Segel system with non-linear diffusion (also referred to as the quasi-linear 
Smoluchowski-Poisson equation) exhibits an interesting threshold phenomenon: 
there is a critical mass M c > such that all the solutions with initial data of 
mass smaller or equal to M c exist globally while the solution blows up in finite 
time for a large class of initial data with mass greater than M c . Unlike in space di- 
mension 2, finite mass self-similar blowing-up solutions arc shown to exist in space 
dimension d > 3. 



1. Introduction 

In space dimension d = 2, the parabolic-elliptic Patlak-Keller-Segel (PKS) system 
is a simplified model which describes the collective motion of cells in the following 
situation: cells diffuse in space and emit a chemical signal, the chemo-attractant, 
which results in the cells attracting each other. If p denotes the density of cells and 
c the concentration of the chemo-attractant, the PKS system reads [13], [19] 

{d t p(t, x) = div [Vp(t, x) - p(t, x) Vc(t, x)} , 
I (t, x) e [0, oo) x R 2 . 

c(t,x) = (E 2 *p)(t,x) , E 2 {x) = -— In [rc| , 

This model may be seen as an elementary brick to understand the aggregation of 
cells in mathematical biology as it exhibits the following interesting and biologically 
relevant feature: there is a critical mass above which the density of cells is expected 
to concentrate near isolated points after a finite time, a property which is related to 
the formation of fruiting bodies in the slime mold Dictyostelium discoideum. Such 
a phenomenon does not take place if the density of cells is too low. More precisely, 
given a non- negative integrable initial condition p with finite second moment, the 
system (CD) has a unique maximal classical solution (p, c) defined on some maximal 
time interval [0, T), T e (0, oo]. Its first component p is non-negative and the mass 
of p (that is, its L 1 -norm) remains constant through time evolution 

||p(t)|| 1= M:=||po||i, te[0,T). 
It is well-known that, if M < 8tt, the solution to ([I]) exists globally in time while 
it blows up in finite time if M > 8ir, see [3] El [HI H2] and the references therein. 
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More recently, it was shown that there is global existence as well for the critical mass 
M = 8 7T, the blowup occurring in infinite time with a profile being a Dirac mass of 
mass 87r pQ. When the mass M is above 8ir, the shape of the finite time blowup is 
not self-similar according to asymptotic expansions computed in [51 [15] (see also [JO] 
for a related problem in a bounded domain). In addition, there is no integrable and 
radially symmetric blowing-up self-similar solution to (TT|) [TBI Theorem 8]. 

In space dimension d > 3, the system (JTJ) seems to be less relevant from the biolog- 
ical point of view as blowup may occur whatever the value of M [SJ [T7] . This means 
that the diffusion is too weak to balance the aggregation resulting from the chemo- 
tactic term. It is however well-known that one can enhance the effect of diffusion to 
prevent crowding by considering a diffusion of porous medium type which increases 
the diffusion of the cells when their density p is large. This is the generalised version 
of the Patlak-Keller-Segel model considered in, e.g., [21 SI 1221 1231 [21]: 



where m > 1, q := l/((d — 2) cr d ), and a d := 2 n d ' 2 /T(d/2) denotes the surface area 
of the sphere S d_1 of M. d . The system also arises in astrophysics [I] (being then 
referred to as the generalised Smoluchowski-Poisson equation), and p and c denote 
the density of particles and the gravitational potential, respectively. 

For ([2]), it turns out that there is only one critical exponent of the non- linear 
diffusion, namely m d := 2(d — l)/d, such that the mass plays a similar role to that 
in flJJ). Indeed, if m > m d the diffusion enhancement is too strong and the solutions 
always exist globally in time whereas if m < m d the diffusion is not strong enough to 
compensate the aggregation term and there are solutions blowing up in finite time 
whatever the value of the mass [221 123]. The relevant diffusion is thus achieved in 
the case when m = m d . In this case, it was proved in [2] that there is a unique 
threshold mass M c > with the following properties: if the mass M = \\po\\i of 
the initial condition p is less or equal to M c , then the corresponding solution to (j5J) 
exists globally in time, whereas given any M > M c there are initial data po with 
mass M such that the corresponding solution blows up in finite time. Thus, for the 
peculiar choice m = m d and d > 3, the system (j2J) exhibits the same qualitative 
behaviour as the PKS system (pQ) in space dimension 2. Still, there is a fundamental 
difference as the latter has no fast-decaying stationary solution with mass 8ir while 
the former has a two-parameter family of non-negative, integrable, and compactly 
supported stationary solutions with mass M c for each d > 3 [2[ Section 3]. 

It is then tempting to figure out whether this striking difference extends above 
the critical mass M c and this leads us to investigate the existence of blowing-up (or 
backward) self-similar solutions with finite mass. More precisely, since mass remains 
unchanged throughout time evolution, we look for solutions (p, c) to (J2j) with m = 
and d > 3 of the form 




d 



(3) 




and 
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for (t,x) G [0,T) x R d with s(t) := [d{T — t)] 1 ^, the time T being an arbitrary 
positive real number. Note that s(t) converges to zero as t increases to the blowup 
time T. 

Our main result is then the following: 

Theorem 1 (Existence of finite mass self-similar blowing-up solutions). There ex- 
ists M 2 G (M c , oo) such that, for any M in (M c , M 2 ], there exists at least a non- 
negative self-similar blowing-up solution (pm,cm) to (j2j) of the form (jHJ) with a 
radially symmetric, compactly supported, and non-increasing profile <3>m satisfying 
\\p M (t)\\i = ||$m||i = M for te [0,T) and \\p M (t)\\oo -> oo ast->T. 

As a consequence of Theorem (TJ we realize that non- negative, integrable, and radi- 
ally symmetric self-similar blowing-up solutions to (j2J) with a non-increasing profile 
only exist below a threshold mass. Another by-product of our analysis is the exis- 
tence of non-negative and non-integrable self-similar blowing-up solutions to see 
Proposition [8] below. 

2. Blowing-up self-similar profiles 

From now on, 

2(d-l) 

d > 3 and m = m d = 



d ' 

and we look for a solution (p, c) to of the form 

with s(t) = [d(T -t)] 1/d and (t,x) G [0,T) x M. d for some given T > 0. We further 
assume that $ enjoys the following properties: 

$ G C(M d ) PI L 1 (IR d ) is radially symmetric and non-negative, 

Inserting the ansatz (j3J) in ([2]) gives that ($, \&) solves 

div (V [$"%)] - $(y)W(y) - $(y) y) = , 
*{y) = (E d *$)(y), 

for y G M d . Since ^ = E^ * $, the radial symmetry of $ ensures that of ^ and, 
introducing the profiles (ip,if)) of ($, \E') 

(5) = ¥>(M) , \P(y) = ^(|y|) , yGM d . 

By [T^l Theorem 9.7, Formula (5)], we have 

1 r i r°° 

for r > 0. We can also write the equation for (p as 

(7) d r (r^ 1 <p(r) d r J(r)) = with J(r) := ^ ~ ^ cp {d - 2)/d (r) - t/j(r) - j 
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for r G (0, oo). Since we are looking for an integrable profile, we formally conclude 
that 

(8) d r J(r) = for r eV v := {s G (0, oo) : <p(s) > 0} . 

In particular, J is constant on any connected component of V v . But, if C is a 
connected component of V<p, we have either 

(9) C = (0, R s ) for some R s G (0, oo] , 
or 

C = (Ri, R s ) for some Ri G (0, oo) and R s G (0, oo] . 

Remark 2. If we additionally assume that the profile tp is non-increasing then 
has only one connected component which is necessarily of the form (jHJ) • 

Now, take a connected component C of V v . It follows from (jHJ) that there is \i G R 
such that 

(10) J(r) = ^ <^~ 2)/d (r) - ^(r) -j = -fi for r G C . 

Owing to the assumed integrability of <E>, the function r h-> r </?(r) belongs to 
L 1 (0,oo) and it follows from (J6j) that the function r i— > r d ~ 2 ip(r) is bounded in C. 
Therefore (jTDJ) only complies with the integrability of $ if R s < oo which implies the 
boundedness of C. Introducing 

and taking the Laplacian of both sides of (ITDJ) yield that H is a positive solution to 

(11) - ^|(r) - — Sr) = fr 2 N (S(r)^ d - 2 )-d) in C, 
with either 

(12) d r E(0) = E(R S ) = if C = (0,i2 s ) 
or 

(13) E(Ri) = E(R S ) = if C = {Ri,R s ). 
A final change of scale, namely 



r/(r):=^H^l , X d := S d ~^ d , & ■= d 1/d 



1/2 



leads us to the following boundary-value problem for rj: either 

+ — ^(r) + r^r)^" 2 ) -1 = 0, r G (0, » d R s ) , 

^(0) = 0, V (ji d R a )=0, 
dr 

(ib) < ^ w + ^ > + " (r) * /M " 1 = ■ r e " A 

T)(n d Ri) = , rj(n d R s ) = . 



(14) 



or 
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We have thus reduced our study to one or several boundary- value problems (depend- 
ing on the number of connected components of V<p) for a nonlinear second order 
differential equation. The purpose of the next section is then a precise study of this 
ordinary differential equation. 

However, before going on, let us point out that (ITT1) is not equivalent to (TTOT) . 
Indeed, since 

d r J(r) = 2(d ~ 1} d r E(r) + f E(s) d ^ 3 d " 1 ds - r , r G C , 

a — 2 r J 

by ©, the fact that H is a solution to (ITT)) only guarantees that d r {r d ~ l d r J(r)) = 
for r G C. Consequently, there are constants C\ and C 2 such that 

^M = -^^, J(r) = £ + reC, 

from which (fTOl) follows only if Ci = 0. On the one hand, if C = (Ri,R s ) with 
< Ri < R s , it is yet unclear whether the boundary conditions (fTBl might imply 
this property. On the other hand, if C = (0, R s ), the boundary conditions (112j) ensure 
that 9 r J(0) = and thus Ci = 0. We shall only deal with this case in the remaining 
of this paper and thus focus on the non-increasing profiles if. 

3. An auxiliary ordinary differential equation 

For a G M, let u(., a) G C 1 ([0, r max (a))) denote the maximal solution to the Cauchy 
problem 

{u"(r, a) + - — - u'(r, a) + \u(r, a) u(r, a) - 1 = , r G [0, r max (a)) , 
r 
ii(0, a) = a , m'(0, a) = , 

with r max (a) G (0, oo] and p = d/(d — 2). 

Clearly, if a = 1 then «(., 1) = 1 is a stationary solution and r max (l) = oo. We first 
show that u(.,a) is global for all a G M and oscillates around the value 1 if a ^ 1. 

Lemma 3. For each a G R \ {1}, r max (a) = oo, and i/ie solution u(.,a) to ( 1161) «s 
an oscillatory function in (0, oo). More precisely, 

• if a > 1, there is an increasing sequence (rj(a))j>o of real numbers such that 
r (a) = 0, 

u'(ri(a),a) = 0, (-l)V(r, a) < /or r G (r^a), r i+1 (o)) , 

u(r 2i (a),a) > u(r 2i+2 (a), a) > 1 > u(r 2i+3 (a) , a) > u(r 2i+1 (a), a) 
for i > 0, 

if a < 1, there is an increasing sequence 0*i(a))t>i of real numbers such that 
n(a) = 

u'(r;(a),a) = 0, (-l) J u'(r, a) < /or r G (r^a), r i+ i(a)) , 

u(r 2i (a),a) > u(r 2i+2 (a) , a) > 1 > u(r 2i+ i(a), a) > w(r 2i -i(a), a) 
for % > 1 . 
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These properties are illustrated in Figure [U Notice that, for a = 7, u(.,7) vanishes 
at a finite r and thus provides a solution to (114j) . 



u(.,a) 




FIGURE 1. Various oscillating behaviours of u(.,a) for a G {0.2, 1,3, 7}. 



Proof of Lemma\^ For any r G [0,r max (a)) consider the functional 
(17) E(r,a) := 1 ' + '_,' n ( r ' a ) 



P 



By (USD, for all r G [0, r max (a)) 



r,a) 



|w'(r,a)| 2 < 



dr * ' r 

Obviously £(r,a) > -p/(p + 1). Owing to (fig) . £(r,a) G [-p/(p + l),£(0,a)] for 
r G [0, r max (a)) which prevents u(.,a) of becoming unbounded at a finite value of 
r, thereby implying that r max (a) = oo. We next argue using Sturm's oscillations 
theorem as in [THl Lemma 9], to establish the oscillatory behaviour of u(.,a) for 
a ^ 1. □ 



According to (114)) . we are interested in finding solutions to the initial value prob- 
lem (fT6l) which are positive and vanish at a finite value of r. We thus focus on the 
case a > and investigate the positivity properties of «(., a). 

Lemma 4. There is a constant a c > 1 swc/i £/ia£ 
• z/ a G (0, a c ) ; i/ien w(r, a) > /or a// r > 0, 
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• if a = a c , then there is R(a c ) > such that 

u(R(a c ),a c ) = 
u'(R(a c ),a c ) = 
u(r, a c ) > forre [0, R(a c )), 

• if a G (a c , oo), then there is R(a) > such that 

u(R(a), a) = 
u'(R(a),a) < 

u(r, a) > for r E [0, R(a)). 
These three possibilities are drawn in Figure EJ 




FIGURE 2. Behaviour of u(., a) for a > a c , a = a c and a < a c . 

Proof of Lemma\J^ For a > 0, we define 

R(a) := inf{i? > : w(r,a)>0 for re{0,R)}. 

Notice that the positivity of a and the continuity of u(., a) guarantee that R(a) > 0. 
We consider the sets 

V := {a > : R{a) = oo} , 

M := {a > : R(a) < oo and u'(R(a), a) < 0} , 
A/" := {a > : #(a) < oo and u'(R(a), a) = 0} . 
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Clearly, V U M U W = (0, oo) and 1 £ V. Actually, if a £ (0, (p + then 
-5(0, a) < and the monotonicity (fl8j) of i? entails that E(r,a) < for all r > 0. 
But, if R(a) < oo, it readily follows from the definition (TlTI) of the functional E 
that E(R(a),a) > whence a contradiction. Therefore, R(a) = oo for any a £ 
(0, (p+ so that 

(19) (0,(J>+1) 1/P )CP. 

Consider now a £ A/"o- Then £7(x) := w(|x|,a) is a radial positive solution to 
the homogeneous Dirichlet-Neumann free boundary problem AU + U p — 1 = in 
B(0,R(a)) with [/ = ^£7 = on dB(0,R(a)). According to [201 Theorem 3 (hi)], 
there is only one value of a for which this solution has a positive radial solution and 
it is unique. Consequently, there is a unique a c > such that J\f = {a c }. 

Consider next a £ A/"UA/o and recall that a > 1 by (fl9"]) . Following [T6l Lemma 11], 
let us assume for contradiction that there is g £ (0,i?(a)) such that u'(g,a) = 0. 
Either g < 1 and we infer from the definition, the monotonicity of E, see ( TT7T) - 
fllSp . and the definition of i?(a) that > E(g,a) > E(R(a),a) > which is a 
contradiction. Or p > 1 and the oscillating behaviour of the solutions implies, using 
the notation of Lemma [31 that g > r 2 (a). This implies that ri(a) < R(a). Then 
u(ri(a), a) £ (0, 1) and using again ffTT|) . f|T8|) . and the definition of -R(a), we conclude 
that > E(ri(a),a) > E(R(a),a) > 0, hence a contradiction. Therefore, 

(20) u'{r,a)<0 for re(0,R{a)) if a £ N U jV . 

Let us now prove that V and M are open subsets of (0, oo) . We first consider a £ M: 
by (|20l there are £> > i?(a) and e > such that u(g,a) < and -u'(r, a) < — 2e for 
r £ (0, £>). By continuous dependence, there is 5 £ (0,a) such that u(g,b) < 
and u'(r,b) < —e for r £ (0, £>) and b £ (a — 5, a + 5). Since w(0,6) = 6 > 0, 
we readily deduce that, for each b £ (a — 5, a + 5), we have R{b) £ (0, £>) with 
u'{R{b), b) < — e < 0. Consequently, (a — 5, a + 5) and A/" is open in (0, oo). Consider 
next a £ V, a > 1. By Lemma [3] and (TTT1) . we have w(r, a) > u(ri(a),a) £ (0,1) 
for r £ [0,r!(a)] and i?(ri(a),a) < 0. By continuous dependence, there is 5 > a 
such that u(r,b) > u(ri(a),a)/2 > for r £ [0,ri(a)], u(ri(a),b) £ (0,1), and 
E(r\(a),b) < for b £ (a — 5, a + 5). Assume now for contradiction that there is 
b £ (a — 5, a + 5) such that R(b) < oo. Owing to ffTTj) . f|T8|) . and the definition of 
i?(6), we obtain > E(ri(a), b) > E(R(b), b) > and a contradiction. Consequently, 
(a — 5, a + 5) CV and V is also open in (0, oo). 

We finally argue as in [TBI Lemma 15] to show that there is A > such that 
(A, oo) C Af. 

Since V and J\f are open subsets of (0,oo), Af = {a c }, (0, (p + l) 1 ^) C V, and 
(A, oo) C A/", we readily conclude that V = (0, a c ) and Af = (a c , oo). □ 



We next study the properties of the map a i— > i?(a). An efficient tool for that 
purpose is the variation of u(., a) with respect to a defined by 

$(r, a) := — (r, a) , (r, a) £ [0, oo) x (0, oo), 
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which solves the second order linear differential equation 

i?"(r, a) + - — - $'(r, a) + p u(r, a) p ~ l #(r, a) = , r G [0, oo) , 

(21) 

■3(0,0) = 1, tf'(0,a) = 0, 
We argue as in [3, [25] to prove the following lemma. 
Lemma 5. If a > a c , there is a unique z(a) G (0,R(a)) such that 

i?(r, a) > /or r G [0, z(a)) , 

0(z(a),a) = 

i?(r, a) < for r e (z(a),R(a)] . 

In addition, u(z(a),a) > 1 and the ratio $(., a)/u(., a) is a decreasing function of r 
on (0,R(a)). 

Proof of Lemma Since the proof follows rather closely that of [25] and [3 Lemma 2.1], 
we sketch it briefly for the sake of completeness. Fix a > a c and set u = u(.,a) and 
$ = #(., a) to simplify notations. We first argue as in [TBI Lemma 17] to show that $ 
vanishes at least once in the interval (0, Zi(a)), where £i(a) denotes the unique zero 
in (0, R(a)) of it — 1. Indeed, ( 1161) also reads 

(u(r) - 1)" + — (u(r) - 1)' + ( u ( r ) - 1) = , r G [0, oo) 

r u{r) — 1 

and (u(r) p — 1)/ (u(r) — 1) < p u(r) p ~ l for r G [0, Zi(a)). It then follows from Sturm's 
comparison theorem that & vanishes at least once in the interval (0, Zi(a)). Let 
z G (0, zi(a)) denote the first zero of "3. 

We now aim at showing that d cannot vanish once more in the interval (z, R(a)). 
To this end, we define 

£( r ) :=r d-i [ u \r) d(r) -u(r) §'(r)} = -r*- 1 u(r) 2 (-} (r) , re[0,R(a)), 



u / 

which encodes the monotonicity of -3/u. It follows from ([IB]) and (121 j) that 

(22) £'(r) = r^ 1 ((p - 1) u p (r) + 1) d(r) , r G [0, fl(a)) . 

Clearly, £'(r) > for r G (0, z) and £(0) = 0, so that £(r) > for r G (0, 2;]. Assume 
now for contradiction that there is G (-2, R{a)) such that 

£(r) > for r G (0, 0) and £(0) = . 

Observing that i?'(z) < 0, we realize that, if #(0) > 0, there is a G (2, g] such 
that i?(r) < for r G (2,0") and $(cr) = 0. In that case, ^X* 7 ) — an d thus 
£(<r) = — cr d_1 u(a) "3' (a) < 0, leading us to a contradiction. Consequently, 

(23) #(g)<0. 
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We next introduce the functions 
2 (u(r)P - 1) 

T r : = 7 — ^ / \p I i Zir) -C{r), 

[p — 1) u{r) p + 1 

((r) := r d [u'(r) 0'(r) + (u(r) p - 1) 0(r)] + (d - 2) r^ 1 u'(r) 0(r) , 
for r G [0,P(a)) and use (05}, (J2U, and ((22) to obtain 
£'(r) = 2 r d " x (w(r) p - 1) 0(r) , 

(24) T'(r) = 2p 2 - -i j «'(r) £(r), 

[(p - 1) u(r)P + 1J 

for r G [0,P(a)). Integrating ( 1241) over (0, and using the negativity of w' and the 
positivity of £ on this interval give 

(25) £(q) = -T(q)>0. 

Since £(g) = 0, we have u(g) &'(g) = u'(g) $(q) and we have 

C(q) = Q(q)M, 

u{g) 

where 

Q(r) := r d [w'(r) 2 + «(r) p+1 - u(r)] + (d - 2) r d ~ l u{r) u\r) , r G [0, R{a)) . 
It then follows from (123]) . fl25l) . and the positivity of u that 

(26) Q(g)<0. 
Finally, define 

P{r) := r d (u'{r) 2 + 2 - 2 «(r)^) + (d - 2) r^ 1 u(r) u'(r) 



p + 1 

for r G [0, R(a)). On the one hand, we notice that 

(27) Pt^Qfrj-titrl-^lt/trf^Qfr), rG[0,P(a)). 
On the other hand, we deduce from (|T6j) and ({TBI) that 

P'(r) = r^ 1 w(r) w(r) p - (d + 2) J , r G [0, P(a)) . 

At this point, we realize that we have necessarily a > (d + 2) (d — l)/(d — 2) and that 
there is s G (0, R(a)) such that P'(r) > if r G (0, s) and P'(r) < if r G (s, P(a)). 
Since P(0) = and P(R(a)) > 0, we conclude that P(g) > and then Q(g) > 
by (|27p . But this contradicts (J26]) . We have thus established that £ is positive in 
(0,R(a)) from which Lemma [5] follows. □ 

We are now in a position to state and prove some properties of the map a i— > R(a). 

Proposition 6. The map a i— »• R(a) is a decreasing function on (a c , oo) and there is 
z\ > such that 

(28) lim R(a) = R(a c ) and lim a {p ~ 1)/2 R(a) = z 1 . 

a\a c a— »oo 
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The monotonicity of a i— ► R(a) is shown in Figure [3] According to numerical 
simulations, the function a i— »• a^ -1 '' 2 R{a) also seems to be a decreasing function of 
a G [a c , oo), see Figured 




FIGURE 3. Monotonicity of the radius i? and a i-> a (p_1)/2 i?(a) (d = 3). 



Proof of Proposition® By Lemma HI u'(R(a),a) < for all a G (a c , oo) and the 
implicit function theorem warrants that R G C 1 ((a c , oo)) with 

dR, , ■d(R(a),a) 
(a) 



<ia u'(R(a),a) 

Since $(i?(a), a) < by Lemma [5j the previous formula implies the strict monotonic- 
ity of a i — > i?(a). We next define 

Ri := sup R(a) G (0, oo] . 

ag(a c , oo) 

If i?/ > R(a c ), there is £ G (R(a c ),Ri) such that u(g,a c ) > by Lemmata [3] and HI 
Then, there is 5 > such that R(a) > ^ for a G (a c , a c + <5). It then follows from the 
continuous dependence of u(.,a) with respect to a and the monotonicity of u(.,a) 
with respect to r that 

= u(R(a c ), a c ) = lim u(R(a c ), a) > lim a) = u(g, a c ) > , 

a\a c a\a c 

and a contradiction. Therefore, i?/ < R(a c ) is finite and we have 

u(R[,a c ) = lim u(R(a),a) = 0, 

a\a c 

from which we conclude that Ri = R(a c ). 
Finally, define 

(29) v(r,a) ■=-""> ( ^1)72 > Q ) » (r, a) G [0, 00) x (0, 00) . 
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Owing to (fTBj) . v(.,a) solves 

v"(r, a) H v'(r, a) + If (r, a)| p_1 v(r, a) — a~ p = , r £ [0, oo) 

r 

u(0,a) = l, u'(0,a) = 0, 
In addition, 

(30) u(r, a) > for r£ (0, a (p " 1)/2 i?(a)) 
for a > a c by Lemma HI Since a _p — > as a — > oo, we have 

(31) lim sup |u(r, a) — u;(r)| = for all g > 0, 



where to denotes the unique solution to 
(32) 



u>"(r) + ~ — - w '( r ) + \ w { r )\ P 1 iu(r) = , r £ [0, oo) 



r 

iw(0) = l, w'(0) = 0. 
By [H], there is z\ > such that 

(33) u>(r) > and w\r) < for r £ [0, zi) , = , it/(zi) < . 

Owing to (1331) . there is 5 > such that w(r) < for r £ (zx, 2! + 5). It then follows 
from fl3~Tj) that, given r £ {z\,Z\ + 5), u(r, a) < for a large enough (depending on 
r), whence a^™ 1 )/ 2 R(a) < r for a large enough by (I3"U1) . Letting r — > Zi guarantees 
that 

limsupa (p ~ 1)/2 R(a) < z 1 . 

Next, if g £ (0,7), we have w(r) > w(g) > for r £ [0,g] and we infer from 
(l3Tj) that u(r, a) > w(g)/2 > for r £ [0, £>] and a large enough. Consequently, 
g < a^" 1 )/ 2 i?(a) for a large enough, from which we conclude that 

liminfa (p - 1)/2 R{a) > z x . 

a— >oo 

Combining the above two inequalities completes the proof of Proposition [61 □ 

The above information allow us to estimate from above and from below a specific 
integral of u(.,a). 

Proposition 7. For a £ [a c , 00), we define 

f R(a) 

M(a) := d \B(0, 1)| / u(r,a) p r d ~ l dr. 
Jo 

Recalling that w is the solution to (1321) and Z\ is its first positive zero, we have 



(34) limM(fl) = M c := d \B(0,1)\ [ 

a^oo J Q 

(35) M 2 := sup M(a) < 00. 

a£[a c ,oo) 



21 



w(r) p r d 1 dr 
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Proof of Proposition^ Let a > a c . Since u(0, a) = a, it follows from the monotonic- 
ity of «(., a) that 

rR(a) 

M{a)<d\B(0,l)\ / a p r d ~ l dr = \B(0, 1)| (a (p " 1)/2 R(a)) . 
Jo 

The upper bound (1351) is then a straightforward consequence of (1281) and the above 
inequality. 

Next, recalling that v (., a) is defined by ( 1291) . we have 

M (a) =d |B(0,1)| / w(r,a) p r d - 1 dr, 

and we infer from (I28p and (l3Tj) that fl34"l) holds true. □ 

4. Proof of Theorem [TJ 

Thanks to the analysis done in the previous sections, we are now in a position to 
construct self-similar blowing-up solutions to ([2]) having either finite or infinite mass. 

Proposition 8. Given a > and T > 0, define 

<p(r) ■= \ d J {d ~ 2) u (fx d r, a) d/[d - 2) for r E [0, oo) if a E (0, a c ) , 

and 

[ \f d - 2) u(/2 d r,a) d /( d -V for r E [0, R(a)/fi d ] 
<f(r) := < if a E [a c , oo) . 

[ for r>R(a)/fj, d , 

Define next ip, and ^ by ^ and (JSJ), respectively. Then the functions (p, c) 
defined by (gj ^n (0,T) x M d mt/i s(t) = [d(T - t)] 1/d a non-negative self-similar 
blowing-up solution to ([2]) wift finite mass if a > a c and infinite mass if a E (0, a c ). 

The proof of Proposition[8]readily follows from the analysis performed in Sections [2] 
and [3j As for Theorem dj it is a straightforward consequence of Proposition [HI the 
threshold values M c and M2 being given by 

M c :=d 1/d ( d _ 2 ) and M 2 := ( rf _ 2 ) M 2 . 

5. Discussion 

We have proved the existence of non-negative, integrable, and radially symmetric 
self-similar blowing-up solutions for . The profile ip of these self-similar solutions 
is compactly supported and non-increasing, and the mass of the corresponding self- 
similar solution ranges in the bounded interval (M C ,M 2 ], the threshold mass M c 
corresponding to the onset of blowup found in [2]. Our analysis thus reveals the 
existence of a second threshold value M 2 > M c of the mass above which no radially 
symmetric and non- increasing self-similar blowing-up solution exist. The meaning 
of this second threshold value for the mass is yet unclear. It is worth mentioning at 
this point that a related situation was uncovered for the critical unstable thin-film 
equation 

d tU = -d x (u n dlu + u n+2 d x u) , (t, x) E [0, 00) x R , 
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in [21] for n £ (0, 3/2). It is likely that, given M £ (M c , Mj\ : there is only a unique 
radially symmetric and non- increasing self-similar blowing-up solution with mass M 
and Figure H] provides some numerical evidence of this fact. Besides this uniqueness 
question, the question of stability of these blowing-up solutions is also of interest. 



Mass 




FIGURE 4. Monotonicity of the mass a i— > M{a). 

Another challenging question is the existence (or non-existence) of integrable pro- 
files tp with a non-connected positivity set as discussed in Section [2j Figure [5] provides 
numerical evidence that, if a > a c is large enough, u(.,a) may have several zeroes and 
each positive "hump" actually corresponds to a solution of (|15|) for suitable values of 
Ri and R s . Whether the additional constraint fTTUl) may be satisfied does not seem 
to be clear. 
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